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Abstract 

A Hamiltonian pair with arbitrary constants is proposed and thus a sort of heredi- 
tary operators is resulted. All the corresponding systems of evolution equations possess 
local bi-Hamiltonian formulation and a special choice of the systems leads to the KdV 
hierarchy. Illustrative examples are given. 

Bi-Hamiltonian formulation is significant for investigating integrable properties of non- 
linear systems of differential equations Q [§] [|| . Many mathematical and physical systems 
have been found to possess such kind of bi-Hamiltonian formulation. There are two im- 
portant problems related to bi-Hamiltonian theory. The one is which kind of systems can 
possess bi-Hamiltonian formulation and the other one is how we construct bi-Hamiltonian 
formulation for a given system if it exists. There has been no complete answer to these 
two problems so far, although a lot of general analysis for bi-Hamiltonian formulation it- 
self has been made. However we can make as many observations on structures of various 
bi-Hamiltonian systems as possible, through which we may eventually find a possible way 
to the final end. 

With such an idea or a motivation to enhance our understanding of bi-Hamiltonian 
formulation, we would like to search for new examples of bi-Hamiltonian systems among 
coupled KdV systems and their higher order partners. There are already some theories 
which allow us to do that. For instance, we can generate soliton hierarchies by using 
decomposable hereditary operators [Q] or by using perturbation around solutions [||. In 
this paper, we would just like to present some new concrete examples to satisfy the Magri 
scheme Q by considering decomposable hereditary operators. 

Let us choose two specific matrix differential operators: 
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<) = jj-, Mj = ad 3 + dj d + 2u ix + imd, u i = u i {x,t) ) < i < N, 



(2) 



1 Email: mawx@cityu.edu. hk 



where ai,Ci,di, < i < TV, are arbitrary constants, but ao / which guarantees the 
invertibility of J. It is known || that J and M constitute a pair of Hamiltonian operators 
with respect to the potential vector u = (uq,ui, ■ ■ ■ ,un) t , that is to say, aJ + bM is a 
Hamiltonian operator for any two constants a and b, which may also be proved directly 
by the Gel'fand-Dorfman algebraic method E| [J7[|. 

Now we can generate a hereditary operator <3? = MJ _1 (see, say, j|| for a general 
proof), since J is invertible. To express this operator explicitly, we need to compute the 
inverse of J. In view of the specific form of J, we can assume 

Pod- 1 p x d~ l ••• PnO- 1 
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where Pi, < i < N, are constants to be determined. It is easy to get that 
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Therefore J J 1 = J 1 J = In+i> where In+i is an identity matrix operator of size (N + 
1) x (iV + 1), leads to an equivalent system of linear algebraic equations for Pi, < i < N: 
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which may be written as 
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Ap = E x , i.e. 






a\ 




Pi 











«o 


a>i 






Pn 










(5) 



The coefficient matrix A is invertible since ao ^ 0, and thus this linear system has a unique 
solution P = A~ l E\. Now we can obtain 



$ = MJ- 1 
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where 

$ f = M,^ 1 = ad 2 + di + 2u ix d~ l + Am, 0<i<N, 
and then the conjugate operator of $ reads as 

p N ^o /?jv_i* + /3iv*i ••• Po^ + Pi^i + --- + Pn^n 
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where 



*i = $J = Ci a 2 + rf. + 2 Ui + 2d~ 1 u i d, 0<i<N. 
Because the Lie derivative of <3> with respect to u x is zero, i.e. 
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we have (see, say, § g g Q) 
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(if m (it + eK n ) - K n (u + eK m )) = 0, K n := $ n u x = 0, m, n > 0. 



e=0 



(10) 

This also implies that a hierarchy of systems of evolution equations ut = K n , n > 0, 
has infinitely many common commuting symmetries {Km}^. All systems in the hierarchy 
have a common recursion operator <3?, since the operator $ is hereditary and has a zero 
Lie derivative with respect to u x : L Ux §> = 0. Moreover, due to the specific forms of < 
i < N, they are all local, although the recursion operator $ is integro-differential. A 
mathematical induction process may easily verify this statement on locality. 

In what follows, we want to show local bi-Hamiltonian formulation for all systems but 
the first one in the hierarchy (note that sometimes systems of soliton equations just have 
one local Hamiltonian formulation in bi-Hamiltonian formulation, such examples can be 
the modified KdV equations and 0(3) chiral field equations [0]). First of all, we observe 
the second system 

u t = Ki = <bu x = MJ~ l u x . 
The vector field J~ 1 u x can be computed as follows 

Pod- 1 p 2 d~ l ••• PnO- 1 
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Evidently we can find or directly prove that 
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where < •, • > denotes the standard inner product ofR w+1 . This means that fo is gradient. 
Now we check whether or not the vector field ^ffo is a gradient field, which is required in 
the Magri scheme O. A direct computation can give 
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, < m < N. 



These equalities yield 
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where the Hamiltonian functional Hi is determined by 

N N 



H X = J Hi dx, Hi=J2 PN-m E A E 



-(CiUjUkxx + diUjUk) + UiUjUk 
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Therefore the system u< = Ki = <fru x has local bi-Hamiltonian formulation 

u t = Ki = <5>u x = J 6 -^ = M 5H ° 
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where Hq and H± are defined by (13) and (15), respectively. 



Secondly, we want to expose bi-Hamiltonian formulation for the other systems ut = 
K n , n > 2. Note that $ = f Ms hereditary, and that /o and ^/o are already gradient. 
According to the Magri scheme M Q, all vector fields ^ n /o, n > 0, are gradient fields, 
namely there exists a hierarchy of functionals H n , n > 0, such that 

fn ■= * n fo = n > 0. (17) 

In fact, the Hamiltonian functionals H n , n > 0, must be equal to 

H n = J H n dx, H n = J < f n (\u),u >dX, n> 0, (18) 
and they are all in involution with respect to either Poisson bracket: 

{H m , H n }j := J -^J-g^ dx = 0, m, n > 0, (19) 
{H m , H n } M := J dx = 0, m, n > 0. (20) 

This way all vector fields K n , n > 1, can be written in two ways as 

K n = $ n u x = $™J/o = .m n f = J-r^, n > 1, 

du 

K n = <b n u x = (J^-'fo = M^-Vq = M— =1, n > 1, 

du 

which provide local bi-Hamiltonian formulation 

ut = K n = $> n u x = J 6 -^ = M^i, n > 1, (21) 
du du 

for all systems ut = K n , n > 1. It follows that the systems ut = K n , n > 0, have common 
commuting symmetries and conserved quantities, which justifies that they constitute a 
typical soliton hierarchy. 

Note that the coefficients appearing in our construction are all arbitrary except the 
requirement of «o 0- Thus the resulting systems may contain many interesting systems. 
A special choice of 

a = c = 1, ati = Ci = 0, 1 < % < N, d { = 0, < i < N, 

leads to the KdV hierarchy under the reduction Ui = 0, 1 < i < N, and thus the above 
resulting systems are called coupled KdV systems. 

Let us now show some examples. Let N = and ao = 1. At this moment, we have 
J = d, M = c d 3 + d d + 2u 0x + Au d, $ = c <9 2 + d + 2u 0x d~ 1 + Au - 
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When Co 7^ 0, do = 0, the corresponding hierarchy is the KdV hierarchy. When cq = do = 0, 
the corresponding hierarchy is a hierarchy of quasi-linear partial differential equations, of 
which the first two nonlinear equations are 



Ut = ®U X 



6u u 0x , u t = ® 2 u x = 30ulu 0x - 



All the vector fields and all the Hamiltonian functionals in this hierarchy are of special 
form cuq 1 uo x and cu™, where c is a constant and m G IN, respectively. 

Let N = 1. The corresponding Hamiltonian pair and hereditary operator become 

c <9 3 + d d + u 0x + 4u d 

cod 3 + dod + uqx + 4uq8 c\d 3 + d\d + u\ x + 4u\d 
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where <5> = c d 2 + d + 2u 0x d 1 + 4u and $1 = ci<9 2 + di + 2ui x 9 1 + 4u\. The first 
nonlinear system is the following 



ut = *>u x = 
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where the Hamiltonian functionals read as 
#0 
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For a general case of N, if we choose 

Qo = 1, aj = 0, 1 < i < N, 
c = l, Ci = 0, < i < iV, 
di = 0, < i < N, 

then the resulting systems are exactly the perturbation systems of the KdV hierarchy 



introduced in Ref. [12] through perturbation around solutions of the KdV equation. 

It should be realized that all first nonlinear systems (ut = K\ = <&u x ) belong to a 
more general class of integrable coupled KdV systems, which was introduced by Giirses 



and Karasu in [13], motivated by the Jordan KdV systems in [14]. Moreover the principle 
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part of our coupled KdV systems, i.e. the systems with di = 0, < i < N , belong to a 
symmetric subclass in the non-degenerate case in [^]. This may be seen by observing the 
coefficients 
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where = 0, i < or i > N, are accepted, after our recursion operators and our coupled 
systems in the case of dj = 0, < i < iV, are rewritten as follows 

N 
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Actually the coefficients defined by (|22| ) satisfy the relations 
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which guarantees [13| that the operators defined by (^) with the coefficients a l - k = 2c* k = 
are recursion operators for the systems determined by ( p4[ ) in the symmetric case of 
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The other nonlinear systems in each hierarchy determined by a hereditary operator $ 
may contain much higher order derivatives of u with respect to x, but they still have a re- 
cursion operator and even bi-Hamiltonian formulation. By taking a scaling transformation 
t — > at, x — > bx, u — > cm, more concrete examples of integrable coupled KdV systems [|l^] 
can be obtained from our systems. 

Compared with the well-known coupled KdV systems (for example, see [15, 16, |T^]), 
the above systems are not really coupled because of the first separated component. The 
Lax pairs or the spectral problems associated with our systems have not been found yet. 
If they are found, master symmetries of the systems can also be presented like ones of the 



well-known coupled KdV systems in [17]. 



Using an idea of extension in Ref. [18], we may obtain much more general Hamiltonian 
pairs starting from the above one. Also we may have other choices, say, 



Ci = consts., 1 < i < 4, 



in Ref. JO]], or start from more general Hamiltonian operator structures to construct new 
integrable systems having bi-Hamiltonian formulation, but we need more techniques in 
manipulation. 
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